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Abstract — This paper presents a theoretical study of oscillatory 
and rotational instabilities of a solid spherical body levitated 
electromagnetically in axisymmetric coils made of coaxial circular 
loops. We apply our previous theory to analyze the static and 
dynamic stability of the sample depending on the AC frequency 
and the position of the sample in the coils for several simple 
configurations. An original analytical approach is introduced 
employing a gauge transformation for the vector potential. First 
we calculate the spring constants which define the frequency of 
small-amplitude oscillations. For static stability the spring con- 
stants must be positive. Dynamic instabilities are characterized by 
critical AC frequencies which, when exceeded, may result either 
in a spin-up or oscillations with increasing amplitude. It is found 
that the critical frequencies increase with the non-uniformity 
of the field. We show that for a spherically harmonic field the 
critical frequency for the spin-up instability in a field of degree I 
coincides with the critical frequency for the oscillatory instability 
in a field of degree / + 1. 



I. Introduction 

ELECTROMAGNETIC levitation melting (ELM) was in- 
vented in the twenties of the last century [1] whereas 
its usage started only at the beginning of the fifties [2] when 
high-frequency power generators became available. The basic 
principle of ELM is simple: a conducting sample, usually 
metallic, weighting from several tens of grams up to several 
kilograms is placed in a coil fed by an AC current with a 
typical frequency ranging from about 10 kHz up to several 
100 kHz. The AC magnetic field induces eddy currents in 
the sample which, in turn, give rise to two effects. On one 
hand, eddy currents interact with those in the coil giving rise 
to a Lorentz force that repels the sample from the coil. On 
the other hand, the Ohmic dissipation due to the induced 
currents provides heating of the sample. In such a way, the 
sample can be levitated and also melted, provided that the 
coil configuration and the current in it are properly adjusted. 
ELM is particularly useful for melting reactive metals with 
high melting points, for example such as Ti, Zr, V, Ta, Mo, 
which often react with the crucible material and so get polluted 
by it [2]. ELM avoids contamination of the melt and allows 
one to carry out soUdification from deeply under-cooled states 

Manuscript received December 13, 2004; revised February 24, 2005. This 
work was supported by Deutsche Forschungsgemeinschaft in frame of the 
Collaborative Research Centre SFB 609 and by the European Commission 
under grant No. GlMA-CT-2002-04046. 

J. Priede was with Forschungszentrum Rossendorf, P.O. Box 510119, 01314 
Dresden, Germany during this work. Presently he is with the Institute of 
Physics, University of Latvia, Miera st. 32, LV-2169 Salaspils, Latvia; 

G. Gerbeth is with Forschungszentrum Rossendoif, MHD Department, P.O. 
Box 510119, 01314 Dresden, Germany 



which is of interest for certain material science applications 
[3]. Besides, ELM is a well-known method for measurements 
of material properties of liquid metals used on ground as well 
as in space [4], [5]. 

The balance of gravity and electromagnetic forces is neces- 
sary but not sufficient for a successful levitation. In addition, 
the sample has to be stable at least to perturbations of suffi- 
ciently small amplitude. First, it means that the reaction force 
due to the displacement of the sample from its equilibrium 
position has to act against that displacement. Otherwise, the 
equilibrium will be statically unstable and the solid sample 
will fall out of the coil or touch it when it is slightly 
perturbed. Similarly, a molten sample may leak out of the 
coil in the result of surface folding and deformation [6]. The 
static stability alone may also be insufficient for a successful 
levitation because sometimes the sample exhibits overstability. 
Namely, there is a restoring force that makes the sample 
to execute oscillations with increasing amplitude so that the 
sample eventually hits the coil or leaves it [2]. Similarly, 
sometimes the sample is observed to spin-up [2], [7], [8]. 
A purely electromagnetic theory trying to account for such 
dynamic instabilities has been proposed in our previous work 
and applied to some simple configurations of magnetic fields 
[9], [10]. In this study, we apply our previous theory to analyze 
both the static and dynamic stability of a spherical sample in 
more realistic axisymmetric magnetic systems made of a set of 
coaxial circular current loops. Although the basic configuration 
is axisymmetric, the perturbation fields are, in general, three- 
dimensional that renders the problem mathematically more 
complicated. To calculate 3D fields for a spherical sample we 
use similar analytic techniques as in Refs. [11], [12]. However, 
our approach differs from the previous ones by an original use 
of the gauge transformation for the 3D vector potential in order 
to satisfy boundary conditions for the induced current. This 
allows us to carry out the analysis only in terms of the vector 
potential without considering the scalar electrostatic potential. 
Therefore our approach is considerably simpler compared to 
[11], [12]. 

The paper is organized as follows. The stability of small 
amplitude oscillations of arbitrary direction is considered in 
Sec. 2 beginning with the governing equations and analytic 
solutions for the vector potential of a circular current loop in 
spherical harmonics. In this section we also derive solutions 
for the axisymmetric base state and the perturbation field 
due to a small displacement. Section 3 presents governing 
equations and analytic solution for the spin-up instability. Nu- 
merical results for both instabilities in several simple inductors 
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are discussed in Section 4 and summarized in Section 5. 

11. Small-amplitude oscillations 

A. Formulation of the problem and governing equations 

Consider a sphere with radius R and conductivity a moving 
at velocity t> in a magnetic field B alternating with circular 
frequency lo. The induced electric field follows from the first 
Maxwell equation as £^ = — V$ — 9* A, where $ is the scalar 
potential of the electric field and A is the vector potential of 
the magnetic field. The density of the electric current induced 
in a moving medium is given by Ohm's law 

j = a{E + VX B) = (t(- V$ - dtA + v xV x A). 

The electric field induced by a translational soUd body motion 
with V being a spatially invariant vector can be represented as 

V xV X A = V{v ■ A) - [v ■ V)A. 

Assuming that the frequency of the alternating magnetic field 
is sufficiently low to neglect the displacement current, the 
second Maxwell equation leads to the following advection- 
diffusion equation 



dtA + (v ■ V)A = V^A 



(1) 



where the gauge invariance of A has been employed to define 
the scalar potential as 

1 



^ = V ■ A 



V A. 



(2) 



Boundary conditions at the surface follow from the continuity 
of the magnetic field 



[Ao\s = [dnAo]g = 



(3) 



where [f]g denotes a jump of the quantity / across the 
boundary S; dn = {n ■ V) is the derivate normal to the 
boundary. 

B. Analytical solution 

1 ) Magnetic field of a circular loop in spherical harmonics: 
The magnetic field of an inductor consisting of a set of coaxial 
circular loops represents a superposition of the fields of the 
separate loops. Therefore, we first consider a single loop lying 
parallel to the x — y plane and centered with respect to the 
z-axis of a Cartesian coordinate system with a sphere at its 
origin, see Fig. 1. Approximation of a toroidal inductor of 
small cross-section by a circular loop is considered in more 
detail in Appendix 1. The loop is suppUed with an AC current 
with amplitude /q alternating harmonically with a circular 
frequency w as /q cos(a;<). Henceforth we use the magnetic 
diffusion time Tm = jJ-ooR^ and the radius of the sphere R as 
time and length scales while the current density and the vector 
potential are scaled by Iq/R^ and iiqIq, respectively. In the 
following, we assume all quantities to be dimensionless with 
the same notation as for the dimensional counterparts used 
so far. Further we consider the vector potential generated by 
the loop A^{r,t) = Ao(r) cos(a)f) with the axisynometric 
ampUtude A^r), where a) = wp.oaR'^ is the dimensionless 




JC V 

Fig. 1. Sketch to the formulation of the problem. 



AC frequency which will be the main parameter throughout 
this paper. 

The amplitude of the current density of such a loop can be 
represented by using the Dirac 5 function: 

3l{r) = e4>5{e - eo)5{r - ro)/ro 

where 6q and ro are the meridional angle and the spherical 
radius of the loop with respect to the center of the sphere 
as shown in Fig. 1. Further, it is advantageous to introduce 
a complex unit vector ^r) — ^ {^x + i^y) that allows us to 
represent the azimuthal unit vector as 

Sfj, = —j= (e^e — Cf^e , 

where and are x and y unit vectors, respectively, and 
the asterisk denotes the complex conjugate. Then the current 
density can be expanded in spherical harmonics as 

Joi'^) = -V2nsm9oS{r - ro)/ro 

OO 1 

xJ^YniM J2 I'^YimiO,^), (4) 

1=1 m=-l 

where Yim{6,^) are spherical harmonics as defined in [14]. 
The expansion coefficients are 



0; m 


?^±1 


ie*; m 


= 1 


ier,; m 


= -1 
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Using the generating function of the spherical harmonics [14] 
the vector potential of a circular loop is obtained as 



> 



■rt 2/ + 1 (r,ro) 

1 



(5) 



where Rq = Tq sin 9q is the cylindrical radius of the loop; 
(r,r')< = min(|r| , (r,r')> = niax(|r| , |r'|). In or- 
der to facilitate the following algebra, it is advantageous 



to introduce the functions X™(r) 



Im \ 



and 



X/"(r) = r''Yim{0,4>), which are the outer and the inner 
solutions of the Laplace equation associated to the spherical 
harmonic Yim{0,(p). Basic properties of those functions are 
given in Appendix 2. For the region occupied by the sphere 
(r < ro) the solution above takes the form 



^21 + 1 

;=1 m=-l 



where 4)q = Q for Tq is assumed. 

2) Axisymmetric basic state: For a sphere at rest the vector 
potential is sought as A(r, = 5? [Ao(r)e*"*] , where is 
complex. For the interior of the sphere, (1) takes the form 



while for the exterior we have 







whereas (2) takes the form of the Coulomb gauge V • Aq 
For the interior, the solution is 

1=1 ^'■^^ 

1 



(6) 



(7) 
= 0. 



(8) 



where gi are unknown coefficients to be determined from the 
boundary conditions; ji{x) is the spherical Bessel function 
of index I [13]. For the exterior, the solution is represented 
as a superposition of external and induced fields Ao(r) = 
Ao(r) + Ao(r) where the latter is sought as 

oo -y-i, s 1 

Alir) = V2nRo ^ -^f^9i E ^'"^/"(r-). (9) 



From the conditions (3) appUed at the surface of the sphere 
at r = 1 we find 



21 + 1 



1 



91 = 



91 



Ji+i\ 



Some useful properties and an efficient algorithm for the 
calculation of gi{x) and its derivative are given in Appendix 
3. 

The time-averaged total force on the sphere is 

Fo = j^{jxB)dV = ^l^?iiijoxB*]dV, 



where the integral is taken over the volume of the sphere V. 
Taking into account that Jo = ^ ^ -^o ^^'^ = Bq + 
Bq, where Bq is a purely real quantity, we obtain Fq = 
h Iv^ [io] xB^dV. According to the momentum conservation 
law the force exerted by a magnetic field on a body is opposite 
to that exerted by the body on the source of the magnetic field 
Fo = — i Jq X 3? [Bq] dV where the integral is taken over 
the space outside the sphere V. The last integral can be taken 
straightforwardly because the external current (4) is defined 
in terms of a Dirac 6 function. Taking into account that the 
induced magnetic field outside the sphere is 

Bi = V2.Ro E E Nrxr 



1=2 

'e;xr 



X / 



m=-2 
■m+1 



X r 



V2Nl^-^ 



V2Nl^+^ 



N, 



i-i 



1/2 



(10) 



where AT™ = ^ ('-'»^)K'+'»)! , we find 

°° / 72 _ 1 

Fo = e,7rE [57-1] //-i/f Ury^T 

where fi = r^'' sin Oq {cos Oq) involves the normaUsed 
Legendre function P™ (x) related to the spherical harmonics as 

loops we have 



Pl^{cos0)e'""l> [13]. In case of N current 



N 

fl = J2 ^-^^n Sin OnPl (cos On) (1 1) 

n=l 

where /„ is the dimensionless current carried by the n-th 
loop; r„, and 0„ are defined analogously as for the single 
loop. Note that the solution (10) alone could be obtained in a 
simpler way as in Ref. [15]. Our following analysis, however, 
is aimed at non-axisymmetric solutions which require a more 
comphcated algebra. 

3} Perturbation due to a small displacement: Consider a 
small displacement x of the center of the sphere from the 
position r to r + x. In the frame of reference related to 
the sphere, this corresponds to a perturbation of the external 
magnetic field 

Al{r + x)^ Al{r) + [x ■ V)A^(r) = A^(r) + \x\ A?(r), 



where 



A?(r) = (e,-V)A^(r) 



1=0 ^^^^ m=-2 

€x = x/ \x\ and 



x/2 



Now the external field represents a superposition of the 
axisymmetric base field AQ{r) and a small but in general 
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non-axisymmetric perturbation A1{r). The perturbation of the 
inner field Ai is governed by the same Eqs. (6, 7) as the base 
field with the difference that now the external field is given 
by Al{r). Consequently, the perturbation of the inner field 
and the corresponding induced field outside the sphere can be 
obtained straightforwardly by replacing the coefficients J™ 
by JJ" in the corresponding solutions (8, 9) for the base field. 
Although the solution of the perturbed field obtained in such 
a way satisfies (6, 7) and the boundary conditions, it turns out 
that the induced field outside does not satisfy the Coulomb 
gauge unless the problem is axisymmetric, i.e., the offset is 
along the symmetry axis. Thus, although the Laplace equation 
is satisfied we formally have V x V x ^4^ = W • A\ = 
j\ ^ which implies that the current is leaking outside the 
sphere. In addition note that the continuity of the tangential 
components of induction following from the continuity of the 
vector potential and its normal derivative at the surface of the 
sphere ensures continuity of the radial, i.e. normal, component 
of the current because rj = r'VxB — 'V - (B x r). Thus, 
the vector potential satisfying the Laplace equation and the 
Coloumb gauge outside the sphere and being continuous at 
the surface ensures that the current is closed in the sphere. 
As shown above, the induced vector potential unambiguously 
follows from the imposed one. Thus, the induced vector 
potential can be modified only by changing the imposed one. 
This can be done by employing the gauge invariance which 
allows us to add to the external vector potential A^r) a 
gradient of some gauge potential defined as 



A^(r) = V2ttRo 



E 

1=1 



21 + 1 



1 



In order to fulfill the Coulomb gauge for the external vector 
potential, Af(r) has to be a harmonic function which is 
reflected in the expression above. So, we obtain a set of free 
coefficients A™ which can be chosen to satisfy the Coulomb 
gauge for the induced field. This gauge transformation results 
in the replacement of the original coefficients J™ by 



e 

V2 



V2 



+ e. 



Then, similarly to the base field, the perturbation of the 
induced field is obtained as 



Aiir) 



\/27ri?o 2^ — —gi 



1=0 



21- 1 



2 

E 

m=-2 



JTXrir). 



From the Coulomb gauge we find 



21 



{21 + 1)1' 



The corresponding induction is 

B\{v) 



1=2 
J 



s/2NV 



m=— 3 

-11 ^ l-l 

V2N^+^ 



NT 



According to our previous theory [10], small -amplitude 
osciUations of the sphere with x{t) = xq cos (fit) give rise 
to a perturbation of the force on the sphere 



Fi = [{-K + inr) ■ xoe'"*] =-K x + r 



dx 
It' 



where K and — J' are the effective electromagnetic stiffness 
and damping tensors, respectively. Further we refer to the 
elements of K and r as the spring constants and growth 
rates, respectively. Note that due to the axial symmetry both 
K and r are purely diagonal matrices with zero azimuthal 
components. At first, the axial synmietry implies that a purely 
axial displacement gives rise to a similar reaction force. 
Second, if the displacement is purely radial the reaction force, 
on the one hand, has to change to the opposite together with 
the displacement because of linearity but, on the other hand, 
it has to rotate around the symmetry axis together with the 
latter Thus, a purely radial displacement causes a purely 
radial reaction force and, therefore, both K and r are purely 
diagonal. Both diagonal elements of K have to be positive 
for the static stability whereas the corresponding growth rates 
have to be negative for dynamic stabiUty. Note that the 
growth rates can become positive and, thus, destabihzing when 
the dimensionless AC frequency exceeds a certain threshold 
depending on the configuration of the field [10]. Both the 
spring constants and the growth rates can be calculated for the 
given direction of the displacement specified by Ex = x/\x\ 
by using the corresponding base and perturbation fields [10]: 



re. 



= -\ JjR[j,]xB'^ + di[jo]xBl]dV 

= j^^[jr]y<BldV. 



Similarly to the integral force, it is advantageous to change 
the region of integration from the sphere to the current loop 



K e 

re 



- = lJ_[foX^[B\]+jlx^[Bi]]dV, 
= -\d^ j_jlx^[B\]dV. 



After some algebra we find 



oo ,2 



(12) 



1=2 
1 



e^xe^x e^(l - y) + e^(ea; • e^) 



Taking into account that 

\ j_3lx BidV = -\ j_3lx {e,-V)BidV 
because j\{r) = (e^ • V)jo(r), where JKr) is defined via 
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Dirac 6 functions we eventually obtain 



21 + 1 



il' - 1) 



^e^ X X e^(l - y) + ez(ex • e^) 



fi-ifi+i fl{l + 2){P-1) 



1/2 



2Z-1 V(2^- l)(2^ + 3) 



Now, we can use this series to calculate both the stiffness 
and the damping coefficients for a given coil defined by the 
coefficients /; according to (11). 

111. Spin-up instability 

A. Governing equations 

We consider a solid sphere, as in Section 2, rotating with 
angular velocity ft in an alternating magnetic field. The 
density of the induced current is given by Ohm's law for a 
moving medium 

j = a{E + v X B) = <T{-dtA - V$ + v x V x A) 

where t; = r2 x r is the velocity of a sohd-body rotation. The 
electric field induced by such a rotation may be represented 
as V xV X A = ^{v ■ A) -{v ■V)A- Axfl. Then from the 
second Maxwell equation with neglected displacement current, 
that corresponds to the quasistationary approximation assumed 
throughout this study, we obtain 



dtA + {v ■ V)A + Axn 



1 



V^A. 



Note that this equation governs the vector potential in the 
laboratory frame of reference where the body rotates while 
the source of the magnetic field is at rest. Further we proceed 
to the frame of reference rotating together with the sphere. In 
a rotating frame of reference, where the sphere is at rest while 
the source of the magnetic field rotates with velocity —v, the 
equation above takes the form as for a body at rest 

dtA{r',t) = —V'A{r',t), 

where r' is a radius vector, which is time-dependent in the 
rotating frame of reference with dtv' = fl x r' . In a rotating 
frame of reference, the body rotation appears as an additional 
time-dependence of the external magnetic field modulating the 
time-dependence of the apphed magnetic field. 
The time-averaged total torque on the sphere is 

According to the conservation of the angular momentum the 
torque exerted by the magnetic field on a body is opposite to 
that exerted by the body on the source of the magnetic field 

M = -i j^rx jgx3?[By dV, 



where the integral is taken over the volume V outside the 
sphere. The last integral can be taken straightforwardly be- 
cause Jo(r') is defined via Dirac S functions. Because (r • 
Jo) = the integral above may be represented as 



M = 



B. Analytical solution 

Further we consider a sphere rotating about an axis per- 
pendicular to the symmetry axis of an axisymmetric external 
magnetic field of a circular loop. The coordinate system related 
to the loop with the z axis as symmetry axis is chosen so that 
the rotation takes place about the x-axis, as shown in Fig. 1. 
It is convenient to analyze this rotation in a coordinate system 
turned about the y-axis by 90 degrees so that the 2;'-axis of the 
transformed coordinate system coincides with the original x- 
axis which is the axis of rotation and, respectively, the x'-axis 
is directed opposite to the z axis of the original coordinate 
system. Then the basis vectors of both coordinate systems are 
related as = — e^, e'y = By, and e'^ — e^- Respectively, 
the complex unity vector is = ^ (e^ ~ ^'rj) + where 
= (e^. + ie'y) . To represent the vector potential of 
a circular loop in the turned coordinate system we use the 
addition theorem for Legendre polynomials [13] written as 



together with the relation 



2l + 3iV0+i 



which yields 

a: 



(r ) - V27ri?o ivj^ V ^ 

X ^ rrxrir') 



m=^—l 



where 



Jm 
I 



1 r^„/ 



NP 



[(e;-e;*)iV;™iF,+i,„(V2,0) 



+ z-^(7V™-iy;+i,„_i(7r/2,0) 
-7V™fy;+i,„,+i(7r/2,0))] . 
Similarly, the current density of the loop can be rewritten as 



jlir') = -\^7rsin6lo 



5{r - ro) 



ro 



1=1 '+1 



Nl I An 



2Z-I-3 



X J2 irYim{e',ci^') 



m=—l 
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In this case, we define a gauge potential 



1=2 



where fi = Rq ^21+°^ -^r- ^J^^\ AJ" are unknown coeffi- 
cients to be determined subsequently from the Coulomb gauge 
for the vector potential induced outside the sphere. Addition 
of VAq to Aq results in the replacement of the coefficients 

jrby 



+e.Nru\\i,)/Nr. 

Now we can simply proceed to the frame of reference rotating 
together with the sphere by redefining the azimuthal angle 
as (j)' = (/>" + (it, where </>" is the azimuthal angle in the 
rotating frame of reference and O is the dimensionless angular 
velocity of rotation. Respectively, the complex unit vector and 
the spherical harmonics change in a rotating frame of reference 
as e'^ = e^'e*"* and X/"(t-') = X/"(T-")e''""*. In a rotating 
frame of reference, the vector potential takes the form 



+ e;'*(e; •X[")e*('"-i)^*) Xj'"(r"). 



Further we consider the current alternating harmonically with 
the dimensionless frequency lo as j'^{r',t) = Jo(r') cos(tDt) 
and the vector potential alternating in the laboratory frame of 
reference, respectively, as 

A'ir', t) = Al{r') cos{u)t) = 3? [Al{r'y^'] . 

In the laboratory frame of reference, the complex ampUtude 
of the vector potential induced outside the sphere is obtained 
as 



,m+l 



1=1 m=-l 

e;*(e;.Xn.9r' 
-e'Ae'^-XT)9T)Xr{r') 



where ,9,™ = ■"+^<A^iZ!) and uj^ 
i(-i(V 

Coulomb gauge V • Aq = we find 

\m T\jm 1 1 



uj + mfi. From the 



ie'* ■ I 



2 (iV™+^] 



V2(iV/"+i)' 



{Nrf _ 

For the field induced outside the sphere we obtain 

00 I 

(r' • V X A^) = inY^lfi E ar^rir') 



x[NJ!l-^Yi:^-\n/2,0) 



The gauge potential Ag vanishes in the expression above and 
so it does not affect the total torque which is 



M 



■■[9l 



X [{l~m + l){l-m + 2)Pll^\0f 
-{l + m + l){l + m + 2)P^l+\0f] . 

Other components of the torque, which are not paral- 
lel to the angular velocity of rotation, vanish because of 
P;"(0)P™=^^(0) = 0. For slow rotations with O «: w we 
have 

C - 5;"" = '^rnndu,9i + O(O^) 
that results in M fa e'CtdcM, where 



I 



(=1 ^ ' ^ ' m=l 



(13) 



X [{l-m+l){l-m + 2)Pl!l^\0f 
-{l + m+ + m + 2)P™+i(0)2] . 

For the rest state of the sphere to be stable to rotational 

perturbations, the torque caused by these perturbations has to 
damp them that corresponds IoOqAI < 0. In the opposite case, 
the rest state is unstable. The state is marginally stable when 
9f2M((Dc) = 0, which is the equation defining the critical 
frequency coc for the spin-up instability. 



IV. Numerical results 

A. Small- ampUtude oscillations 

In the following we present some numerical results for 
an inductor consisting of two circular loops of equal radii 
placed coaxially at the distance H from the mid-plane. The 
current having the same amphtude may flow either in the same 
or in opposite directions in both loops. Subsequently, both 
these cases are referred to as symmetric and antisymmetric, 
respectively. We start with a simple symmetric case when 
both loops coincide {H = 0) making up a single loop with 
doubled current amplitude. The forces along with axial and 
radial spring constants versus the distance of the sphere from 
such a loop of dimensionless radius i? = 2 are shown in 
Fig. 2 for various dimensionless frequencies. As seen in Fig. 
2(a), at large enough distances the force, which is repulsing, 
raises with decreasing distance between the loop and the 
sphere. At some distance, the force attains a maximum and 
tends to zero as the sphere approaches the center of the loop. 
Such a reducing force implies that the position of the sphere 
is statically unstable to axial perturbations. The position is 
statically stable to axial perturbations when the axial spring 
constants shown in Fig. 2{b) are positive that requires the 
sphere to be placed further away from the loop than the point 
of maximum force which is at z ~ 0.75. A position being 
axially stable may become unstable radially when the distance 
from the loop becomes too large (see Fig. 2b). Thus the range 
of statically stable positions for a single loop is rather limited 
(0.75 < z < 1.5) and does not change significantly with the 
frequency of the magnetic field. 
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Further let us turn briefly to the static stabiUty of a sphere 
in an arrangement of two coaxial current loops of the same 
radius placed symmetrically at distance H from the center of 
the sphere. For the sake of brevity we restrict the static stabiUty 
analysis to the perfect-conductor approximation a) = oo. Axial 
and radial spring constants are shown in Fig. 3 versus the 
distance H at various radii of the loops for both symmetric and 
antisymmetric arrangements. In the symmetric arrangement, 
which is usually used for heating, the position of the sphere is 
statically unstable in radial direction when the loops are too 
much separated, and it becomes unstable in the axial direction 
when the loops are moved too tightly together. Overlapping of 
the stability ranges for radial and axial perturbations depends 
on the radius of the loops. For i? = 2 the range of axial 
stability ends at iJ w 1 where the range of radial stability 
begins. Thus, there is no overlapping of the stabiUty ranges 
in this case that impUes a static instability regardless of the 



distance between the loops. A range of positions staticaUy 
stable to both axial and radial perturbations is possible only for 
sufficiently small radu of the loops. For instance at R = 1.25 
this range is approximately 0.5 < 2: < 0.75. In contrast to 
the symmetric arrangement, an antisymmetric one, which is 
usually used for positioning of the sample, ensures a statically 
stable state to radial perturbations regardless of the distance 
between the loops (see Fig. 3b on right). The position of the 
sphere is stable to axial perturbations for any radius of the 
loop provided the loops are separated by a distance larger 
than the radius of the sphere (H > 0.5). There might be an 
axial instability for smaller separations of the loops at < 2 
(see Fig. 3a on right). 

Forces together with axial and radial spring constants versus 
the distance of the sphere from the mid-plane of both symmet- 
ric and antisymmetric arrangements of two loops are shown in 
Fig. 4 for various distances between the loops in the perfect- 
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Fig. 3. Axial (a) and radial (b) .spring constants versus the hialf-distance between the loops for both symmetric (on left) and antisymmetric (on right) 
arrangements at various radii of the loops in perfect-conductor approximation (to = oo). 



conductor approximation (u ~ oo). It is seen that the position 
of the sphere in the middle -point (z = 0) of a symmetric 
arrangement (Fig. 4 on left) is unstable either in radial or axial 
directions for almost all distances between the coils. However, 
there is one distance between the coils, « 1, at which 
both spring constants cross zero almost simultaneously. Thus, 
this distance is an optimal one for the static stability of the 
sphere in the middle-point. Note that the distance H = R/2 
provides the most uniform magnetic field in the vicinity of the 
middle-point of the symmetric arrangement. This is because 
the next to the leading order contribution for the magnetic 
field of symmetric arrangement vanishes at the aspect ratio 
H/R = cot 6*0 = 0.5, i.e., PaHcos^o) = in (5). Similarly, 
the distance H = providing the most uniform gradient 

of the magnetic field in the vicinity of the middle-point of an 
antisymmetric arrangement is expected to ensure the widest 
range of statically stable positions for loops of large enough 
radius. 

The dynamic stability of small-amplitude radial and axial 
oscillations is characterized by the corresponding growth rates 
Frr and F^z defined by (13). The threshold of linear stability 
is given by the frequency above which the corresponding 
growth rate becomes positive. At the threshold frequency, the 
amplitude of oscillations turns from decaying to growing in 



time. A characteristic feature of this instability is that the 
growth rate, becoming positive above the threshold frequency, 
reaches a maximum at a certain frequency and tends to zero 
at higher frequencies. The frequency at which the growth rate 
attains a maximum is referred to as the most dangerous one. 
Threshold frequencies along with the most dangerous one and 
the corresponding maximal growth rates are plotted in Fig. 
5 versus the half-distance between the loops at various radii 
of the loops for both axial and radial oscillations of a sphere 
about the middle-point of both symmetric and antisymmetric 
arrangements of two loops. The most dangerous frequency is 
always above the threshold one. As seen in Figs. 5(a) and 5(^), 
at large enough distances between the loops, the threshold and 
the most dangerous frequencies tend to Wc = 27.682 and ivd — 
47.196 in the synnmetric arrangement, and to lUc = 11.609 
and u)d = 18.792 in the antisymmetric one, respectively. Note 
that these critical frequencies for the oscillatory instability 
in spatially quadratic and Unear magnetic fields coincide 
with those for the spin-up instability in linear and uniform 
magnetic fields, respectively, which will be considered later. 
For each radius of a symmetric arrangement, there is a certain 
optimal distance between the loops that maximizes the critical 
frequency. For large enough radii the optimal distance between 
the loops is approximately equal to the radius, 2H = R, 
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(on right) arrangements of two coaxial loops of radius R = 2 and various distances from the mid-plane at a> = oo. 
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which is optimal also for the static stability as considered 
above. With decreasing of the loop radius, the optimal distance 
becomes shghtly larger (see Figs. 5a and 5b on left). In the 
antisymmetric case, the critical frequency changes weakly on 
reducing of the distance H until the loops approach a distance 
H w 0.75 to the mid-plane (see Figs. 5a and 5b on right). 
Further approaching of the loops results in the increase of the 
critical frequency which is the larger, the smaller the radius. 
This raise is caused by the nonuniformily of the magnetic field 
in the vicinity of the loop. Although the critical frequency 
for axial oscillations is, in general, shghtly lower than that 
for the radial oscillations, there is no significant difference 
between both. The maximal growth rate for the symmetric case 
changes weakly with reducing distance between the loops until 
the distance becomes comparable to the radius of the sphere 
(z w 0.5) (see Figs. 5c and 5d on left). Closer approaching of 
the loops causes a significant increase of the growth rate for 
both axial and radial oscillations. Again, the smaller the radius 
of the loops, the larger the increase of the growth rate. In the 
antisymmetric case, the growth rate reaches a maximum as 
the loops approach to the mid-plane at a distance comparable 
to the radius of the sphere (z « 1) (see Figs. 5c and 5d on 
right). 

B. Spin-up instability 

Spin-up instability occurs when the frequency of the alter- 
nating magnetic field exceeds the threshold ujc defined above. 
The spin-up rate being proportional to OqM becomes positive 
as the frequency raises over the threshold and attains a max- 
imum at some higher frequency ujd subsequently referred to 
as the most dangerous one. These critical frequencies together 
with the corresponding maximal spin-up rates for a sphere 
at the middle-point between symmetric and antisymmetric 
arrangements of two loops are plotted in Fig. 6 versus the half- 
distance between the loops of various radii. As seen in Fig. 
6(a), at large enough distances between the loops the threshold 
and the most dangerous frequencies tend to Wc = 11.609 
and uici = 18.792 in the symmetric arrangement, and to 
ujc = 27.682 and ujii = 47.196 in the antisymmetric one, 
which correspond to the spatially uniform and hnear fields, 
respectively. In the symmetric case, the critical frequencies 
have a minimum at the distance between the loops equal 
to their radius: 2H = R (see Fig. 6a on left). Note that 
this is converse to the critical frequencies for the oscillatory 
instability considered previously which have a maximum at 
this distance. This minimum is due to the maximal uniformity 
of the field in the vicinity of the middle-point achieved at this 
distance between the loops. Critical frequencies increase with 
the non-uniformity of the field which becomes particularly 
significant when the loops are approached closer than the 
radius of the sphere (H » 0.5). In this case, the smaller the 
radius of the loops, the higher the critical frequencies. The 
spin-up rate raises as the loops are approached and attains a 
maximum at i7 w 0.5. For the loops with large enough radius 
the maximum of the spin-up rate is attained when both loops 
merge together forming a single one (see Fig. 6b on left). In 
the antisymmetric case, the critical frequencies begin to raise 



significantly when the loops approach to the mid-plane at a 
distance of the radius of the sphere: « 1 (see Fig. 6a on 
right). As seen in Fig. 6b on right, the spin-up rate attains a 
maximum approximately at the same distance. 

V. Summary and conclusions 

This work presents an analytic stability analysis of a solid 
sphere levitated electromagnetically in an axisymmetric mag- 
netic field induced by a set of coaxial circular loops. Stability 
is considered with respect to both oscillations of small amph- 
tude and arbitrary direction and rotations perpendicular to the 
axis of symmetry of the field. Both oscillations and rotations 
are found to become growing in time when the dimensionless 
frequency of the field exceeds a certain critical threshold u)c 
depending on the configuration of the field for each kind 
of instability. The growth rates of both instabilities become 
positive as the frequency is raised over the corresponding 
threshold, attain a maximum at some higher frequency and 
tend to zero at high frequencies. The frequency at which 
the maximum is attained, is referred to as the most dangerous 
one for the corresponding instability. Critical frequencies and 
the corresponding maximal growth rates are calculated for 
arrangements of two loops of the same radius carrying the 
same total current which may flow either in the same or in 
opposite directions. These arrangements, corresponding to the 
usual heating and positioning fields, are referred to as the 
symmetric and antisymmetric ones. It is found that critical 
frequencies for both oscillatory and rotational instabihties 
increase with the nonuniformity of the magnetic field. The 
lowest dimensionless critical frequencies are u)c = 11.609 
and LUd = 18.792 which are the same for both the spin-up 
instability in a uniform magnetic field and the oscillatory in- 
stability in a spatially linear field. This coincidence represents 
a particular case of a more general rule: the critical frequencies 
for the spin-up instabiUty in a spherically harmonic field of 
degree I coincide with the critical frequencies for oscillatory 
instability in a spherically harmonic field of degree / -|- 1. 
Thus, the critical frequencies Uc = 27.682 and uid = 47.196 
for the spin-up instability in a linear field coincide with the 
critical frequencies for the oscillatory instability in a spatially 
quadratic magnetic field which is the case at the mid-point of 
a symmetric arrangement of two loops except the case that 
the distance between the loops is equal to the radius of the 
loops. In this particular case, the quadratic term vanishes and 
the instability is dominated by the fourth-order term which has 
critical frequencies Wc = 85.252 and = 158.6. Thus, the 
distance 2H = R is optimal in order to avoid the oscillatory 
instability of the sphere at the mid-point of the symmetric 
arrangement. This distance, which ensures the most uniform 
magnetic field at the middle-point, yields a minimum of the 
critical frequency for the spin-up instability. Note that there 
is no similar optimal distance for the oscillatory instabihty in 
antisymmetric arrangement where the critical frequency can 
be increased only by the nonuniformity of the magnetic field 
when the loops are approached close to the surface of the 
sphere. Hence, the stabilization of a sphere by an optimal 
design of the inductor is rather hmited and active means of 
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Fig. 6. Critical (Wc) and most dangerous (uJd) dimensionless AC frequencies (a) and maximal growth rates (h) for a spin-up instability in symmetric (on 
left) and antisymmetric (on right) arrangements of two loops versus the half-distance of the loops with various radii. 



Stabilization, like an additional steady magnetic field, may be 
necessary. 

Note that the maximal dimensionless growth rates for both 
oscillatory and spin-up instabilities are small, typically F ~ 
10~^. From the physical point of view we have introduced 
the growth rate as a viscous-type friction coefficient at the 
velocity. Thus the physical dimension of F is s/m and, 
consequently, its scale is Tm/R, because we have used the 
magnetic diffusion time and the radius of the sample R 
as time and length scales, respectively. In order to obtain the 
actual friction force we have to take the product of velocity, 
friction coefficient and the magnitude of the characteristic total 
electromagnetic force Fq which usually is comparable to the 
gravity of the sample mg, where m is the mass of the sample 
and g the free fall acceleration. By comparing the negative 
effective friction force to the inertia we obtain an estimate of 
the characteristic growth time of the instability tq ~ p^^- 
which for a characteristic size R ^ 10~^ m and a conductivity 
a ~ 10^ f2~^m~^ of the sample leads to tq ~ 10'* s. Thus, 
although the magnetic field has a dynamically destabilizing 
effect, the development of the instability is expected to be 
very slow and, thus, hardly observable. In conclusion, note 
that the frequency of the magnetic field is the only parameter 
determining the threshold of dynamic instabilities as long as 



no external damping, for instance due to a surrounding gas, 
is taken into account. In the opposite case, the threshold of 
instability would depend not only on the frequency but also 
on the amplitude of the current. 

In conclusion note that the dynamic instabilities resulting 
from the effect of motion of a conducting body in an AC 
magnetic field can be interpreted from an alternative physical 
point of view which can explain the weakness of this insta- 
bility and also suggests other possible instability mechanisms. 
In the given AC magnetic field, the electromagnetic force on 
the spherical body at rest depends only on its position. If the 
position changes, it takes some time for the electromagnetic 
force to relax to its time-averaged stationary value for the 
given position. This delay is caused by the finite time of 
magnetic diffusion in the body which, taking place over the 
skin depth, is thus comparable to the AC oscillation period t. 
The electromagnetic force at instant t on the body in motion 
may be represented as F(r{t — r)). Assuming a sufficiently 
high AC frequency, we can expand this representation in a 
power series of small r yielding F{r{t — r)) « F{r{t) — 
Tv) F(r{t)) — t{v ■ V)F. Thus, the delay of the magnetic 
field results in a force component proportional to the velocity 
which is obviously analogous to a viscous friction force. Since 
for a statically stable position the effective electromagnetic 
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reaction force {x ■ 'V)F is directed against the displacement 
X, the effective friction force —t{v ■ V)l^ is directed with the 
motion. Consequently, the delay of the magnetic field results 
in an effective electromagnetic friction force with a negative 
and, thus, destabiUzing friction coefficient. It is important to 
note that the friction coefficient is small because it involves 
the delay time which is short and comparable to the AC period 
when the relaxation of the magnetic field is solely due to 
the electromagnetic diffusion over the skin depth. The same 
arguments may be extended to the coupling of the magnetic 
field and the temperature distribution in the body caused 
by its temperature-dependent electrical conductivity. In this 
case, the delay of the magnetic field would be dominated by 
the thermal relaxation time which is much longer than the 
electromagnetic diffusion time and, thus, expected to cause 
stronger dynamic instabilities than the pure electromagnetic 
mechanism considered here. In addition, the decrease of the 
electrical conductivity for increasing temperature, which is 
caused by the Joule heating of the AC magnetic field, and the 
related reduction of the non-dimensional frequency u) might 
be important in levitation experiments. 

Appendix I 

Approximation of a toroidal inductor of small 
cross-section by a circular loop 

Consider a toroidal inductor of arbitrary cross-section car- 
rying an axisymmetric and purely azimuthal current. The 
magnetic field of such an inductor may be represented as a 
superposition of the fields of separate circular loops constitut- 
ing the inductor, and the corresponding vector potential may 
be written as the integral over the inductor cross-section S: 
A{r) = Jg Ao{r,r')ds' , where 

Ao{r,r) = — / 

An J„ \r-r'\ 

is the nondimensionalized vector potential of a single circular 
loop defined by the radius vector r' = r'ej, + z'e'^ with r' 
and z' being the radius and axial position of the loop in 
cylindrical coordinates. Subsequently assume the skin effect to 
be negligible as if the inductor would be fed by a direct current 
driven by the gradient of an electrostatic potential $. An 
axisymmetric and purely azimuthal current distribution with 
density j = je^ implies $ to depend solely on the azimuthal 
angle </>. Consequently, j = -V$ = e^^^, where || = C 
is a constant that can be related to the total current in the 
inductor /q = effective radius 

Js r 

as C = loR/ S. Further we assume the inductor cross-section 
to be located about some position tq = tqC + zqC^, which 
will be specified later, and approximate the field distribution 
by the following multipole type expansion: 



Airy 



,IoR 
'AttS 



I [1 + (r' - ro) • Vo] 
Js 



ds'dcj)' 



where the operator Vq acts on rg. Now, if we define ro as the 
position of the mass center of the cross-section ro = ^ rds, 
the second dipole-like term cancels in the expression above 
which, thus, reduces to 



A{r) 



Iro e> 
47r Jo \r - 



2- e'^d^ 



(15) 



The last expression defines the vector potential of a cir- 
cular loop located at ro and carrying the effective current 
/ = loR/ro- For example, in the case of an inductor 
represented by a torus with major and minor radii ro and 
ri, the effective radius can easily be found from (14) as 



i? = f 



1 + 



while the mass center of the 



cross-section coincides with its geometrical center. Note that 
the solution (15) is accurate up to the quadrapole-like term 
neglected in (15) the magnitude of which relative to the re- 
maining term may be estimated as ~ |r' — ro|^ / |r — ro|^ ~ 
(d/i?)2, where d is the characteristic size of the cross-section 
and R is the characteristic distance to the cross-section center. 
Obviously, the neglected term becomes significant only at 
distances comparable to d. 

Appendix II 

Properties of the functions xp{r) and X™(r). 

Similarly to Ref. [11] we also use a complex unity vec- 
tor = {ex + iey) , where and By are the corre- 
sponding Cartesian unit vectors, and employ the outer and 
inner solutions of the scalar Laplace equation, Xp{r) = 
r"'"^yi,m(6',<A) and X/"(r) = r%rn{0,(l)) associated with 
the spherical harmonic Yi,„,(6.(t)) [13]. These variables pro- 
vide a simple algebra for calculation of the gradient operator 
V = J-;ir + e*|- -I- e^lj used throughout this study: 



dxi 
dxi 

dz 
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1 Nr 



drj 
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2Z- 17V™/ 



-X 



where the asterisk denotes the complex conjugate and NJ^ = 

^ ^'""i+i"'"'"^' - Th^ corresponding relations for the com- 
plex conjugate functions follow straightforwardly from the 
above relations and the property of spherical harmonics: 
Yi*m{^A) = (-l)"'>l,-m(6',(/)). Another relation used for 
the calculation of the torque in Sec. 3 is: 



r X VXp = i 



{l-m)X, 



m+1 



■ ro\ 



.^/2 A^, 

e^NT 
V2 N," 



+ e^mXT 
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-{l + m)X- 
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Note that our definition of Xp is slightly different from that 
used in [11] leading to a bit more complicated algebra but 
simpler resulting expressions. 

Appendix III 
Calculation of the coefficients gi{z). 

The coefficients gi{z) and g'i{z) in (10,13, 13) can efficiently 
be calculated using continued fractions [16]. For this purpose 
we rewrite gi{z) = = hi{z)hi^i{z) where hi{z) = 

•'V^^^^ and ji{z) is the spherical Bessel function of index I. 
Further, applying a recurrence relation for the spherical Bessel 
function of index l+l we obtain the continued fraction: 



' {2l + 3)/z-hi+i{z) 
_ 1 1 1 

~ {21 + 3)/z- {21 + 5)/z- {21 + 1)1 z- ■■" 

allowing us to calculate hi{z) provided that is 
known. But asymplolic properties of Bessel functions suggest 
that for large index I hi{z) ^ ^. Thus for any z we can 
choose sufficiently large I' and then truncate the fraction by 
approximating hi' (z) by the previous expression that allows us 
to calculate back the necessary hi{z). Note that the forward 
recurrence for hi{z) is not practically apphcable because it 
is numerically unstable similarly to its counterpart for Bessel 
functions. The other necessary quantity g[{z) = contained 
in (13,13), can be calculated in a similar way by expressing it 
as g'i{z) = {21 + l){hf_j^{z) - gi{z))/z. For a) » 1 we have 
gi{yju}/i) ^ —1 + 2^^ which follows from the corresponding 
asymptotics of the Bessel functions [13]. 
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